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1. Introduction 


As found by Griineisen about fifty years ago the thermal expansion is at low 
temperatures proportional to the energy content and deviates at high tempera- 
tures from this proportionality towards higher expansion. For these phenomena 
there has for a very long time been a rather loose connection between the results 
of theoretical attempts and those of experiments. The present writer has made 
a re-examination of these and related phenomena on the basis of a phenomeno- 
logical analysis of the experimental results combined with partly new assumptions 
concerning the sources of energy and expansion (1). In a recent paper (2) the 
energy content of solid Au, Ag, Cu, and Al was analyzed in detail. This analysis 
will now be followed by a discussion of the thermal expansion of the same metals 
including also the expansion at the melting point and in the liquid state. 


2. The energy content 


We shall first briefly recapitulate the results of the recent investigation (2) 
concerning the energy content of the metals Au, Ag, Cu, and Al in their solid 
state. It was shown that the total observable energy content per mole £,4,, obtained 
by integrating over the specific heat from absolute zero to a given temperature, 
can be described as a sum of four parts 


Eops = Ey + Ha + £1 + Ex. (1) 


E,» is the observable part of the energy due to vibrations of the atoms. The 
key to the dividing of the energy is the assumption that this observable energy 
of vibrations at high temperatures approaches the value 
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where ZH, is the vibrational zero-point energy, which is not directly observable. 
This means that the sum of zero-point energy and thermal vibrational energy 
approaches the classical value 3 #7’, R being the gas constant. At low tempera- 
tures equation (2) gives too small values for Ey», which should reach zero at 
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T=E,/3 R. The difference between the real Hy, and that calculated from equa- 
tion (2) is a quantity which was called the residual zero-point energy Ey, and 
which with increasing 7’ passes from EH, to 0. The full expression for E, thus is 
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HE. in equation (1) is the thermal energy of the free electrons. 
E, and Ey could be described as functions of the temperature by the expres- 
sions ; 
E,= A; exp (— H;/RT) (4) 
and 
Ey = Ay exp (— Hy/RT) (5) 


where A,, H;, Ay and Hy are empirical constants. These expressions reveal the 
sources of the energies H; and Hy to be thermal lattice defects. The order of 
magnitude of the heat of formation Hy is 1 eV and the corresponding defects are 
identified as vacancies. The order of H, is 0.2 eV and the structure of these 
defects is still not established. 


3. Assumptions concerning the thermal expansion of solid metals 


As was first found by Griineisen there is at low and medium temperatures a 
good proportionality between the increase of volume and the increase of energy 
with increasing temperature. At these temperatures the energy is predominatingly 
of vibrational origin and the thermal expansion is therefore proportional to Ey. 
We now assume that the vibrational part of the expansion remains proportional 
to EH, also at high temperatures. This proportionality also follows from conside- 
rations of the thermal density fluctuations (1). At high temperatures H,,, is de- 
scribed by equation (2), which means that the vibrational part of the expansion 
varies linearly with 7. 

We have no basis for a judging if the small electron energy H,, which is 
proportional to 7”, has any influence on the thermal expansion. We shall try the 
assumption that an eventual influence is negligible. 

The structural energy Z; most probably originates from some lattice defects and 
it is natural to assume a contribution to the expansion that is proportional to 
Hy, as long as the number of defects is small. 

The energy Hy most probably originates from vacancies with a heat of for- 
mation Hy. Their number pro mole is Ey/Hy. If each vacancy increases the 
volume in the same way as an added atom, the relative increase of volume due 
to vacancies should also be Hy/Hy. According to calculations by Eshelby (3) and 
Jongenburger (4) the real volume expansion should be smaller by a factor of 0.6, 
owing to relaxation of the lattice. As, however, these calculations are based 
on a static lattice, it is not obvious whether they also apply to the thermally 
agitated and defective lattice which may be expected at the temperatures where 
vacancies appear at an appreciable extent. We therefore prefer first to try the 
simple assumption that the relative increase of volume is equal to Hy/HAy. 


Summing up we arrive at a tentative expression for the relative increase of 
volume which may be written 
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For practical reasons the difference between the volume V at a given tempe- 
rature and the volume V, at absolute zero is here taken relative to the volume 
V, at 0°C, which is the reference temperature for most dilatometric measure- 
ments. This choice of V, instead of V, as reference value has an influence on 
the constants a and b of the order of one per cent. The last term in equation 
(6) gives the fraction of atoms replaced by vacancies. It is of course not rational 
to compare this fraction with the change of volume referred to the volume Vy at 


90°C. On the other hand it is difficult to say what should be most rational, a 


reference to V, or to the volume expanded by the lattice vibrations. This makes 
an uncertainty of about 5% which, however, does not matter as the present 
experimental uncertainty about the term due to vacancies will be found to be 
much larger than that. 

For a comparison of equation (6) with experimental results concerning (V—V.)/ 
V, for Au, Ag, Cu, and Al we take Ey, E1, Ey, and Hy from the results of the 
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1023 6130 5524 
1073 6524 5813 
1123 6929 6102 
1173 7346 6391 


earlier paper (2). The constants a and 6 thus are the only unknown quantities 
to be determined. The vibrational term with the constant a dominates at tempe- 
ratures below and around room temperature. The term due to the primary defects 
with the constant b becomes important at elevated temperatures. The last term 
due to the vacancies appears at still higher temperatures and has no very great 
influence on the determination of a and 6. 


4, Comparison with experimental results 


The experimental and calculated values for Au, Ag, Cu, and Al needed for a 
test of equation (6) are put together in the Tables 1 to 4. The volume expansions 
are calculated from the dilatometric measurements by Nix and MacNair (5) on 
Au, Cu, and Al in the range from —180 to 700, 500 and 400°C, respectively, by 
Esser and Eusterbrock (6) on Au, Ag, Cu, and Al from room temperature to the 
neighbourhood of the melting points and by Austin (7) on Au in the same range. 
Nix and MacNair have published extensive tables of direct observations from 
which we have obtained the values in the tables by interpolation and transfor- 
mation from linear to volume expansion. Esser and Eusterbrock have published 
smoothed values for every fiftieth degree and Austin for every hundredth degree. 

Ey is for high temperatures calculated by equation (2) with the values of E, 
taken from the previous paper and introduced in the headings of the tables. For 


low temperatures the values of Ey, are taken directly from the tables of the 
previous paper. 
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E, and Ey are calculated by means of equations (4) and (5), with the empirical 
values of A;, H;, Ay, and Hy from the earlier paper (2) shown in the headings of 
the present tables. 

The empirical constants a and 6 are determined by trial and successive approxi- 
mations. The values of V, are obtained from extrapolations to L,;,=0. The last 
columns of the tables show the differences between the mean values of (V — V,)/Vo 
from the different series of measurements and the sum of the three terms on the 
right side of equation (6). The periodicity of these errors, especially marked in 
the table for Ag, are partly due to the fact that the expansion values are taken 
from results smoothed by the use of power series of 7’. 

Close to the melting points the observed values are definitely smaller than the 
calculated ones. This discrepancy cannot be satisfactorily removed by introducing, 
according to Eshelby (3), a factor smaller than unity in the vacancy term of 
equation (6). There is perhaps a possibility that the specimens have yielded to 
the slight pressure of the measuring arrangements and it is useless to discuss 
details concerning the ranges below the melting points until the question is made 
quite clear on the experimental side. 

The proportionality constants a, b, and | /Hy of equation (6) for the relative 
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Table 4. Thermal expansion of solid Al. 


Experimental data from Nix and MacNair (N. M.) and Esser and Eusterbrock (E. E.). 
H,=730, A, = 4040, Hy= 4720, A, = 89000, H; = 13200 cal/mole. 
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Table 5. Empirical proportionality constants to equation (6). 


10°xa 10°xb 10°/H,, 
Au 0.702 1.65 4.3 
Ag 0.97 1.60 3.6 
Cu 0.865 1.82 5.3 
Al 1.21 2.40 7.6 


expansions connected with the vibrational energy, the energy of the primary 
defects, and the energy of the vacancies, respectively, are put together in Table 
5. These constants are of the same order of magnitude but still markedly differ- 
ent. An interesting observation is that b is about the double of a. Now the first 
term in equation (6) could also be written 2a U,,, where Uy, is the potential half 
of #\,. This means that the proportionality factor between thermal expansion 
and potential energy is about the same for vibrations and primary defects. 


5. The thermal expansion illustrating the separate contributions 
from vibrations and defects 


For the range of temperature below and not too high above room temperature 
we can in equation (6) omit the last term due to the vacancies and introduce 
the value of HE, from equation (3). Starting from 0°C with 7 =T,, V=Vo 
E., = Ez, and Ey= Ey, we obtain the equation 
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Fig. 1. Separation of contributions to the thermal expansion from vibrations and defects. Data 
calculated from the results of Nix and MacNair. 
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—a3 R(T —T;) =4 (Be, — Bayo) +b (Ey — Ey). (7) 


Nix and MacNair have published detailed series of observations covering the range 
of temperature where the residual zero-point energy decreases to zero and the 
structural energy HH; increases from zero. Values of the expression on the left 
side of equation (7) derived from their measurements on Au, Cu, and Al and the 
values of the empirical constant a from table 5, are introduced in figure 1. For 
Al and Cu the two terms on the right side of equation (7) overlap. For Au the 
two contributions are well separated. There is a range of temperature on both 
sides of room temperature where the expansion coefficient of Au is practically 
constant and equal to a 3 R. 


Table 6. Ratios of the changes of volume and energy on melting compared % 
corresponding ratios for the primary defects. 


WH cealimpley ws. eo ORS | 3060 | 2690 | 3110 2500 
10° A Vn/Vo A Em, (cal/mole)-? . 1.76 1.50 1.42 2.52 
10° A V/V A Ey (cal/mole)-?. . . 1.65 1.60 1.82 2.40 


6. Nature of the expansion at the melting point 


As shown in table 5 the expansions due to vibrations, primary defects, and 
vacancies are characterized by differing values of the ratio AV/V, AE of rela- 
tive expansion and increase of energy. It is therefore of fundamental interest for 
an understanding of the melting phenomenon to compare this ratio for the changes 
at the melting point with those for the other phenomena. The data for a com- 
parison are put together in table 6. The values of the relative change of volume 
A Vin/Vms, Where Vs is the volume of the solid at the melting point, are all 
taken from a critical compilation by Kubaschewski (8) where especially the A V,,- 
values were chosen between rather differing experimental values. The heat 
of fusion is from the compilation by Kelley (9). The ratios obtained come very 
near the ratios b for the primary defects obtained from table 5, and introduced 
in the last row of table 6. The greatest difference, in the case of Cu, cannot at 
present be judged as a real discrepancy outside the possible uncertainty of the 
experiments. We may conclude that the melting of the metals in question is 
mainly a discontinuous increase of the primary defects. 

Another interesting question that now presents itself is to what extent the 
liquid state is saturated as to primary defects. The energy of formation of pri- 
mary defects gives the change of heat when one mole of atoms passes over from 
the undisturbed lattice to the defects and is obtained from measurements in a 
range of temperature where the concentration of defects is still small. It seems 
probable that H; will decrease at higher concentrations and that the observed 
H, therefore constitutes an upper limit for the heat of saturation. The ratio 
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where £,,, is the energy due to primary defects just below the melting point, 
and AE,, the heat of fusion, therefore means a lower limit for the degree of 
saturation. From the data given above we obtain for Au, Ag, Cu, and Al x =0.77, 
0.66, 0.72, and 0.60, respectively. We conclude that these metals are rather 


highly (perhaps fully) saturated with respect to the primary defects in the liquid 
state. 
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al ‘str CAl/mole degree 
| 108 me (cal/mole)-? + eo 4.3 5.0 4.1 8.9 
[oti (caljmele)~27. 37. 2... 4.3 3.6 5.3 7.6 


7. Nature of the expansion in the liquid state 
Also for the liquid state it is, in principle, possible to obtain information about 


4 the nature of the thermal expansion from a study of the ratio between the struc- 
tural part of the expansion coefficient 


- _1av_ 

bs Oster = V, a7 a3RkT (9) 
3 and the structural part of the atomic heat capacity 

x 


Cuw=C—-3R-y Tn (10) 


- where C is the observed molar heat capacity and y 7’, the heat capacity of the 

_ free electrons at the melting point. Though the accuracy of the determinations 
of dV/dT and C for the liquid metals are still very unsatisfactory we have for 
a rough orientation collected the available data in table 7. The relative changes 

of volume per degree dV/Vmed T, where Vy, is the volume of the liquid at the 

melting point, are taken from the compilation by Kubaschewski (8) and the heat 
capacity in the liquid state from the compilation by Kelley (9). The values used 

for y are the same as in the previous paper (2). The last line is taken from the 
last column of table 5. 

From the last two lines of the table we see that the ratio between the changes 
of volume and heat with temperature is of the same order for the liquid state 
as for the formation of vacancies. In fact the uncertainty of the compiled data 
is such that the ratios may well be identical. We conclude that the structural 
part of the thermal expansion of these liquid metals is mainly due to the for- 
mation of vacancies. This is in good agreement with our earlier conclusions that 
the liquids are already rather saturated as to primary defects. The concentration 
of vacancies remains low high above the melting point. 


Stockholm, Royal Institute of Technology. 
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